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Abstract. Permanent magnets with high-permeability yokes have been widely used 
in watt balances for supplying a robust and strong magnetic field at the coil position. 
Subjected to the mechanical realization, only several Br{z) (radial magnetic field along 
the vertical direction) profiles can be measured by coils for field characterization. In 
this article, we present an algorithm that can construct the global magnetic field of 
the air gap based on N{N > 1) additional measurements of Br{z) profiles. The 
proposed algorithm is realized by polynomially estimating the Bz(r) function with 
analysis of basic relations between two magnetic components in air gap, i.e., Br(r^ z) 
and Bz{r,z), following the Maxwell’s equations. The algorithm, verified by FEM 
simulations, can characterize the three-dimensional contribution of the magnetic field 
for a watt balance magnet with acceptable accuracy, which would supply basic field 
parameters for alignment and misalignment corrections. 
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1. Introduction 

The watt balance, which was originally proposed by Dr. B P Kibble at the National 
Physical Laboratory (NPL, UK) in 1975 [1], is an experiment currently for precision 
measurement of the Planck constant h, and in future for maintaining one of the seven 
base units, the kilogram [2], The operation of a watt balance is divided into two 
independent measurement modes, respectively known as the weighing mode and the 
velocity mode. In the weighing mode, a magnetic force generated by current-carrying 
coil in a magnetic held is balanced by the weight of a test mass m as mg = BLI, 
where g denotes the local gravitational acceleration, B the magnetic hux density at coil 
position, L the wire length of the coil, and I the current in coil. In the velocity mode, 
the coil is moved in the magnetic held with a velocity v in vertical direction, inducing a 
voltage U = BLv on the coil terminals. By a combination of the weighing and velocity 
modes, the geometrical factor BL is eliminated and a virtual watt balance equation, i.e., 
UI = mgv, would be obtained. As a result, the Planck constant h can be related to the 
test mass m by comparing the electrical power to the mechanical power in conjunction 
with the quantum Hall ehect [3] and the Josephson ehect [4]. Since the proposal in 1975, 
the watt balance has been widely spread and persuaded at many national metrology 
institutes (NMIs) [5-13]. The detailed principle and recent progress of watt balance 
experiments at NMIs can be found in several review papers, e.g., [14-16]. 

In order to obtain enough magnetic force, e.g., 5N, and in the meanwhile keep the 
power assumption of the coil as low as possible, a strong magnetic held is preferred in 
watt balance. In realization, the watt balance usually employs permanent magnets as 
the magnetic hux source and uses soft yokes to guide the generated hux through a small 
air gap where the coil is placed. As most of the magneto motive force (MMF) in such 
magnetic circuits would drop across the air gap, the magnetic hux in the air gap would be 
strong. Besides, the held boundary on yoke-air surfaces is sharp with high-permeability 
yokes, and hence the held uniformity in the air gap is good. But one disadvantage for 
such magnetic circuit is that the magnetic held gradient in the horizontal direction is 
obviously increased, which may bring difficulties for the alignment [17]. As is noticed, 
the alignment is one of the most difficult procedures in operating a watt balance, and 
several techniques that can relax the alignment have been proposed: the simultaneous 
measurement of two operating modes has been carried out in watt balance at the Bureau 
International des Poids et Mesures (BIPM) [8], which has been further developed by Dr. 
I A Robinson in [18]; new experimental designs of watt balances, of course, can relax 
the alignment [19]. 

Here in this paper, a misalignment correction idea based on measurement of coil 
position and representation of magnetic held in the air gap, is considered. The merit 
of the correction method is that it would relax the alignment and can be applied on 
all the existing watt balances without redeveloping new apparatus. The coil position in 
both modes of a watt balance can be well measured as a necessary requirement of the 
experiment. In order to eliminate the misalignment error to a level below 1 x 10“® by 
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Figure 1. A typical magnetic structure used in watt balances (sectional view). Two 
circles in the upper and lower parts of the magnet denote the main flux in the magnetic 
circuit. The flux direction shown is one of two opposite cases when the northern poles 
of permanent magnets are facing the inner yoke. 


correction, a key procedure is to know the three-dimensional (3D) magnetic prohle in the 
operating interval of the air gap. However, limited to the mechanical realization, only 
several vertical prohles of Br{z) at different horizontal coordinates can be determined 
by the gradient coil (GC) method [20] or a high resolution magnetic probe. Knowing a 
set of subjected information, a natural question is: how can we solve the whole hied in 
the air gap by using minimum measurable quantities? To answer this question, in this 
paper, we present an algorithm that can represent the 3D magnetic prohle of the air 
gap based on N{N > 1) additional measurements of Br{z) prohles. 

The outline for the rest of this paper is organized as follows: the algorithm principle 
is presented in section 2 which is originated from polynomial estimation of the Bz{r) 
function following Maxwell’s equations in the static magnetic circuit, and in section 3, 
numeral verihcations by a cubic estimation of the Bz{r) function with one additional 
Br{z) measurement are shown; in section 4, the equations for misalignment correction 
are reviewed and discussed; a conclusion is drawn in section 5. 

2. Principle 

2.1. Magnet structure 

A typical magnetic structure employed in watt balance with one coil and two permanent 
magnet rings has been shown in hgure 1. The magnetic circuit was hrst presented in the 
BIPM watt balance group [21], and later followed by the METAS Mark II watt balance 
at the Federal Institute of Metrology, Switzerland [22], the NIST-4 watt balance at 
the National Institute of Standards and Technology, USA [20], the MSL watt balance 
at the Measurement Standards Laboratory, New Zealand [23], and the KRISS watt 
balance at the Korea Research Institute of Standards and Science, South Korea [13]. 
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The shown magnet strncture is preferred due to the advantages including generation of 
a strong magnetic field in the air gap, a good magnetic shielding, and the possibility of 
a compact size realization. 

In the shown symmetrical magnetic structure, two permanent magnet rings are 
arranged with opposite magnetization poles and their flux, guided by high-permeability 
yokes, runs horizontally through the air gap where the coil is suspended. The air gap 
is located between two parallel cylindrical yokes, i.e., the inner yoke and the outer 
yoke, and its width is conventionally several centimeters. The height of the air gap is 
typically tens of centimeters and the central part with several centimeters is the applied 
measurement interval in the velocity mode. Normally, both the weighing and velocity 
modes are operated in a chosen measurement interval along the central Bj.{z) prohle, 
which can be well adjusted with a relative hied change of several parts in 10^. 

Note that there are some other realizations of watt balance magnets, e.g., [11,24]. 
Although the model and analysis presented in this paper is based on the magnet shown 
in hgure 1, they can be also applied in other magnet structures. 


2.2. Magnetostatic equations 

We hrst present the generalized magnetostatic equations to describe the held in the air 
gap for further discussion of the proposed algorithm. For a static magnet circuit shown 
in hgure 1, the diherential forms of Maxwell’s equations for the magnetic held in the air 
gap are written respectively as 


V • B = 0, 

(1) 

V X B = poJ, 

(2) 


where J denotes a three dimensional vector of the current density, B is the vector of 
the magnetic hux density in the air gap, and p-o is the permeability in vacuum; V is the 
Del operator dehned as V = ro{d/dr) + 0o(l/r)((9/90) -|- zo{d/dz) where (ro, 0o, zq) 
is a unit vector in the cylindrical coordinate (r, (p, z). 

Writing the vector B in equation (1) in forms of three components in the cylindrical 
coordinate, i.e., B = {Br, B^p, Bz), we obtain 

l^rBr{r,(p,z)] ^ ldBp,{r,(p,z) ^ dBz{r,(p,z) 
r 


+ 


= 0 . 


( 3 ) 


dr r d(j) dz 

Since the presented magnet in hgure 1 has a. r — z symmetrical structure, it meets the 
following relations for the magnetic hux density in the air gap that 


dBr{r,(j),z) dBz{r,(j),z) 

= ^ = ^ = 

As a result, equation (3) can be rewritten in a two dimensional form as 

ld[rBr{r, z)] ^ dBz{r,z) 
r dr dz 


= 0 . 


( 4 ) 


( 5 ) 
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Similarly, if J = ( J^j,, Jz), equation (2) can be expressed by three components in 
the cylindrical coordinate respectively as 


ldBz{r,((),z) 
r dcj) 


dB^{r,(t),z) 

dz 


/ioJr(r,0,2r), 


( 6 ) 


dBr{r,(j),z) 

dz 


dBz{r,(l),z) 

dr 


/ioJ0(r,0,2r), 


l d[rB^{r,(j),z)] _ 1 dBrjr, (j), z) 
r dr r dcj) 


(7) 

( 8 ) 


Note that equations (6)-(8) are written in the space of the air gap, where the current 
density should equal zero, i.e., Jr = J<i> = Jz = 0. Besides, based on equation (4), both 
left sides of equation (6) and equation (8) equal zero, therefore, equations (6)-(8) can 
be simplihed by an equation in two dimensional form as 
dBr{r,z) _ dBz{r,z) _ 

dz dr ' ^ ’ 

Equations (5) and (9), which present the basic relations between two components 
of the magnetic flux density, i.e., Br{r, z) and i?^(r, z), are the magnetostatic equations 
for the magnetic held in the air gap of the shown watt balance magnet. 


2.3. Fringe field effect of the air gap 

In this section, some discussions on the fringe held ehect of the air gap are prepared for 
the algorithm. Based on the electromagnetic knowledge, the vertical component of the 
magnetic hux density in the air gap, i.e., Bz{r, z), is generated by a hnite height of the 
air gap, which is known as the fringe hied ehect. The schematic of the fringe held for 
the magnetic hied in the air gap has been shown in hgure 2(a), and it is known that 
the fringe ehect would bend the magnetic hux at both upper and lower ends of the air 
gap. As a result, the absolute value for the vertical magnetic component Bz{r,z) will 
increase along a vertical axis departing the center position {z = 0). In theory, if the 
upper and lower magnets are symmetrical, only the magnetic held line at the radical 
center has zero vertical magnetic component, i.e., Bz{r,0) = 0. 

As high-permeability yokes are used, the boundary on the yoke-air surface is very 
sharp, and it is reasonable to consider the magnetic hux at both yoke-air surfaces to be 
ideally radical, i.e., Bz{rc — 6q/2^z) = Bzir^ -|- (5o/2, z) = 0 where 6o is the air gap width. 
It can be proved in mathematics, as well as in conjunction with Maxwell’s equations, the 
Bz{r) function is a typical one-extremum function. A typical Bz{r) function curve with 
z ^ 0 has been shown in hgure 2(b), which usually has two approximate zero values at 
surfaces and one extremum Bz{rm, z) around the central air gap. 

2.4- Algorithm 

The proposed algorithm employs A^-|-1 measurements of Br{z) prohles at diherent radii, 
i.e., Br{ri, z) where z = 1, 2,..., A^, A^ -|- 1. The central prohle Br{rc, z) (rc is the central 
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Figure 2. (a) Schematic of the fringe field effect for the magnetic filed in the air gap. 
(b) A typical (r) profile along the radical direction between two surfaces of the inner 
and outer yokes at a certainty vertical position z{z ^ 0). 

radius of the air gap), conventionally being measured in the velocity mode of a watt 
balance, is included and indexed as i = iV + 1. The other iV measurements of Br{z) 
prohles indexed form 1 to iV are specially designed for supplying additional information 
for the algorithm. It should be noted that a larger number of N will obtain a better 
estimation accuracy but it would obviously increase the complexity of the measurements, 
e.g., a multi-coil with N upper coil and N lower coils is required if the GC method is 
applied. The goal of this paper is to represent the whole air gap magnetic held by very 
small numbers of additional Br{z) measurements, e.g., iV = 1, with acceptable accuracy 
for relaxing the alignment. 

In the calculation, iV -|- 1 measurements of Br{z) prohles are taken, and hence 
N + 1 values of dB^^r, z)/dr along the radical direction of the air gap can be calculated 
by equation (9). As discussed in section 2.3, it is reasonable to make the vertical 
component of the magnetic hux density at both yoke-air surfaces equal zero, i.e., 
Bz{tc — 5o/2, z) = Bz{rc + 6o/2,z) = 0. When N is a. small number, the information 
with knowing -|- 1 diherential values and 2 initial values at yoke-air surfaces for the 
Bz{r) function, however, is still not enough for solving the magnetic hux density along 
the radical direction. 

In order to make up the information lacking and ensure the held representation 
accuracy, a pre-estimation of Bz{r) is applied. It has been pointed out in section 2.3 
and shown in hgure 2(b) that the Bz{r) function is a one-extremum function with two 
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zero values at yoke-air surfaces. A polynomial estimation of Bz{r) function should be 
feasible and can be applied in the algorithm. The polynomial estimation employs the 
Br{rc,z) prohle and the other N additional Br{z) measurements at different radii for 
best estimating the vertical magnetic component in the radical direction, which can 
yield a polynomial estimator with maximum order K = N + 2 . The main idea of the 
proposed approach is to £t the Bz{r) function at every vertical coordinate by a -|- 2 
order polynomial function in dimension r of the air gap as 

r ^ N 

z) = a{z){r - -7 y)(r - - j) JJ(r - /3i), (10) 

i=l 

where a(z) is a gain factor; (3i is the (z-|-2)th solution of B^^r, z) = 0 while /3_i = rc-|-(5o/2 
and /So = — 6o/2 are considered as its hrst and second solutions. 

It can be seen from equation (10) that Bz{r,z) is expressed by a product of two 
functions with separated directions, i.e., the two dimensions (r and z) of the Bz{r,z) 
function is decoupled, therefore, the calculation of the magnetic held in the air gap is 
greatly simplihed. For each B^{ri, z) measurement, it can be written by equation (9) as 


dBr{ri,z) 


dB^{ri,z) 


i = 1,2, ...,N,c. 


( 11 ) 


dz dr 

The left part of equation (11) is directly measured while the right part of equation 
(11) can be expressed by the estimator following equation (10). As a result, A^ -|- 1 
unknown parameters, i.e., f3i{i = 1, 2,..., N) and a, can be solved by A^ -|- 1 equations in 
(11). Substituting the solved parameters back into equation (10), the Bz{r,z) function 
on the whole air gap is obtained. 

Knowing the Bz{r, z) function in the air gap, we can directly calculate values of the 
differential function dBz{r, z)/dz at any single point {r,z). According to equation (5), 
d[rBr{r, z)]/dr is solved as 

d[rBr{r,z)] _ dB^{r,z) 

dr dz ' ^ ’ 

To ensure the calculation accuracy in the actual operation area of a watt balance, 
the initial condition of rcBr{rc, z) is applied in the algorithm, the rB^{r, z) function and 
hence the Br{r, z) function in the air gap is solved. 

In summary, both Bz{r,z) and Br{r,z), i.e., a 3D magnetic held representation in 
the air gap, can be solved by polynomial estimation of Bz{r) function. The minimum 
number of additional Br{z) measurements for the estimation equation (10) in theory 
can be zero, i.e., A^ = 0. Note that this case is established only when a B^{z) prohle 
with r 7^ Tc is known. For watt balance experiments, the Br{rc, z) prohle is preferred to 
obtain a larger vertical interval with good uniformity but the gain factor a in this case 
cannot be solved. Therefore, in order to represent the 3D magnetic held in the air gap 
for watt balance, at least one additional Br{z) prohle measurement should be taken. 
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Figure 3. Parameters setup of the numeral simulation (unit:mm). 


3. Numeral verification 

3.1. Simulation setup 

In order to evaluating the model accuracy presented in section 2, here some numeral 
simulations based on hnite element method (FEM) are performed. In the simulation, 
the geometrical parameters of the magnet have been shown in hgure 3. The relative 
permeability of the yoke is set as = 1000 and the magnetic strength of the permanent 
magnet is set as 800kA/m in the vertical .2 direction. Note that in the simulation of a 
watt balance magnet, the nonlinearity of the yoke permeability, instead of a constant 
number, should be considered. However, in our approach, the nonlinear information 
has already been contained in measurements of Br{ri,z) prohles. In such cases, the real 
magnet can be simplihed by a constant-permeability one and no obvious error would 
be generated in evaluating the model accuracy. In the simulation. Eve magnetic prohles 
Hj.(219mm, 2 ;), Br{215mm, z), Br{210mm,z), Br{205mm,z), and i?j.(201mm, 2 ) shown 
in hgure 4, are hrstly calculated as possible known conditions of the presented model to 
simulated the actual measurements of either a GC coil or a magnetic probe. 

3.2. A case with iV = 1 

Here the simplest condition with iV = 1 additional Br{z) prohle when the third-order 
polynomial estimator is calculated as an example. In this case, Bz{r,z) is considered 
to be a cubic function of r with arbitrary 2 ; position. Following equation (10), Bz{r,z) 
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Figure 4. Br{z) profiles applied in the calculation of the presented model. 


function can be expressed as 




B,{r, z) = a{z){r - r, + —)(r - - —)(r - /3i), 


and its partial derivative is written as 
dB^ir, z) 


"Oi 


dr 


= 2(r - rc)(r - I3i)a{z) + [(r - Vc) - 


Using the profile Br{rc,z), a(z) can be solved as 

4 dBr{rc,z) 


a{z) = 


Si 


dz 


(13) 


(14) 


(15) 


Here an additional measurement i?r(^"c + S,z) where 6 G {—5q/2,5q/2) is taken. 
Based on equation (14), it can be obtained that 


dBr{rc + S, z) 
dz 


A2 

= 5{2rc + 3(5 — 2j3i)a{z) — -^a^z). 


Substituting equation (15) into equation (16), /5i can be calculated as 

_ M ^ _ dBr{rc + S,z)/dz 

^ 2 86 dBr(ri.,z)ldz 


(16) 


(17) 


In the calculation, only two Br{z) prohles, Br{rc, z) and Br{ri, z) shown in hgure 4, 
are applied. An area with r G (195mm, 225mm), G (—50mm, 50mm), where the watt 
balance is operated, has been focused. The calculation results of Br{r,z) and Bz{r,z) 
with Tc = 210mm and ri = 201mm have been shown in hgure 5. 

It can be seen from the calculation results that the amplitude of Bj.{r) decays along 
r direction approximately following a 1/r relation. This relationship can be physically 
explained: the total hux through the air gap, <h, is a hxed number only determined by 
the MMF of the permanent magnet and the magnetic reluctance of the circuit; since the 
magnetic hux on each surface is approximate uniform, the total hux can be expressed as 
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Figure 5. Calculation results of magnetic flux density distribution with vertical 
profiles Br{210mm, z) and _Br(201mm, z). (a) Magnetic flux density distribution of 
Br{r,z) in 3D. (b) Magnetic flux density distribution of Bz{r,z) in 3D. (c) Top view 
of the magnetic flux density distribution Br{r,z). (d) Top view of the magnetic flux 
density distribution Bz{r, z). 


$ = 27rrdB{r) where d is the height of the air gap. Therefore, Br{r) = ^/{27rrd), i.e., 
Br{r) decays along r direction with 1/r. Besides, the calculation of Bz{z) also clearly 
shows the fringe hied effect, i.e., B^ increases quickly in the air gap when 2 ; is departing 
from the central radical surface z = 0. 

In order to show the calculation accuracy of the proposed algorithm, the errors in 
the focused space, compared to FEM simulations, are calculated. The error maps for 
both Br{r, z) and B^{r, z) with different choices of the additional z) prohle, i.e., ri 

respectively equals to 201mm, 205mm, 215mm, and 219mm, have been shown in hgure 
6. It can be seen from hgure 6 that the calculation error is closely related to the choice 
of the additional Br{ri,z) prohle. The maximum relative error in diherent cases when 
= 1 is about 0.1%-0.4% for Br{r, z) and about 4%-16% for Bz{r, z). The calculation 
result shows that the error when the additional B^{z) is taken at ri > Tc is smaller than 
that when the additional Br{z) is taken at ri < Tc. Therefore, it is suggested to take 
an additional B^{z) measurement in the outer half of the air yoke, e.g., r = + 2(5o/3, 

when a cubic estimator is applied. 
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Figure 6. Filed calculation error in Gauss when N = 1. The additional Br{z) profile 
is adopted at the radius ri that (a) ri = 201mm, (b) ri = 205mm, (c) ri = 215mm, 
and (d) ri = 219mm. In each subgraph, the left is the magnetic flux density error of 
Br{r, z) while the right is the magnetic flux density error of Bz{r, z). 
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3.3. High order estimators 

Another simulation example when = 2 is also taken and the simulation results have 
been shown in hgure 7. Two valuable conclusions have been found from the calculation: 
Firstly, two additional measurements, Br{ri, z) and Br{r 2 , z), can not be set symmetrical 
on Tc, i.e., ri — Tc ^ r^ — r 2 . Otherwise, an infinite value for the gain factor a would 
obtained. In this case, a is expressed as 

^ _ ^2^^ d[Brir.,z)-Br[r,,z)] ^ 

where 5i = ri— rc and 62 = r 2 — rc. Obviously, the condition (5i + ^2 = 0 is not allowed. 
Secondly, the calculation error with ri and r 2 chosen in the same half side of the air gap, 
either outer or inner parts, is smaller than the cases with ri and r 2 in different halves of 
the air gap. Therefore, if the measurement with iV = 2 is taken, it is better to choose a 
same half of the air gap in radical direction, e.g., ri = rc — Sq/S and r 2 = — 2(5o/3, or 

ri = rc + 5o/3 and r 2 = rc + 2(5o/3. 

Higher order estimators with N{N > 2) measurements of Br{z) prohles contain 
more hied information, and in theory can improve the calculation accuracy. However, 
it would synchronously increase the complexity for both measurement and analytical 
solving unknown parameters of the estimator. Limited to the paper length, no further 
calculation with higher order estimators is presented. 


4. Discussion 


The misalignment error ^ of a watt balance has been presented in [17], expressed as 
UI 


g = 


1 = 




l ^y'^y | 


^ x^x 






(19) 

FzVz FzVz ' F^Vz ' FzVz ' FzVz ' FzVz ’ 

where F = {Fx;, Fy, Fz) and T = (ra;,rj^,r 2 ) are magnetic forces and torques generated 
in the weighing mode. Note that F and T are results of the interaction between the coil 
current I and the magnetic hux '0, which can be written respectively as 


r - -TC^ j- 


( 20 ) 

( 21 ) 


In above equations, 6 = {9x,0y,9z) is the rotation angle of coil, and (i,j,k) is a unit 
vector in (x, y, z) space. 

Base on measurements of Br{z) prohles and the presented model, 3D functions 
of the magnetic hux in the magnet air gap ^|J = {gpx^i’y^i’z) can be calculated, and 
hence the magnetic forces F = {F^, Fy, Fz) and torques T = (ra,,ry,r^) as functions 
of geometrical space are known. Further, in conjunction with optical measurements of 
the velocity v = {vx,Vy,Vz) and the angular velocity to = (t0x,t0y,t0z), the misalignment 
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Figure 7. Filed calculation error in Gauss when N = 2. Additional Br{z) profiles 
are adopted at two radii, ri and r 2 . (a) ri = 201inm, r 2 = 205mm; (b) ri = 215mm, 
ri = 219mm; (c) ri = 205mm, r 2 = 219mm; (d) ri = 215mm, r 2 = 201mm. In each 
subgraph, the left is the magnetic flux density error of Br{r,z) while the right is the 
magnetic flux density error of Bz{r, z). 
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error in equation (19) can be determined, and a correction of misalignment errors in 
theory can be applied. 

It should be noted that all terms on the right side of equation (19) are tinny values 
compared to 1, which conventionally should be adjusted below 1 x 10“® by alignment. As 
presented, the proposed algorithm is expected to represent the 3D magnetic prohle with 
an accuracy of several parts in 10^ with few additional measurements of Bj.{z) prohles. 
On the one hand, under such situations, a whole held parameter can be supplied for 
alignment analysis. On the other hand, if the held information is used for misalignment 
correction, the alignment would be relaxed, e.g., several parts in 10^. 

5. Conclusion 

An algorithm employing a polynomial estimator for functional ht of Bz{r) is presented 
to characterize the 3D magnetic held of a watt balance. Based on Maxwell’s equations, 
the relations of two magnetic components in air gap of a watt balance magnet, Bj.{r, z) 
and Bz{r,z), are directly modelled by decoupling the Bz{r,z) function in r and 2 ; 
dimensions. The accuracy of the presented model is evaluated by comparisons to the 
FEM simulation, and the calculation results showed the alignment of a watt balance 
experiment can be greatly relaxed by few additional Bj.{z) prohles. As the model relies 
on real measurements of Br{z) functions, the nonlinear information of the yokes has 
been contained in the calculation, and in theory can be applied in watt balances with 
diherent magnetic prohles. The model supplies basic magnetic held parameters, which 
can be used for alignment and misalignment corrections. 
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